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摘 要 : 多 层 网 络 的 同步 性 研究 是 网 络 科 学 研究 中 一 个 重要 的 课题 , 其 中 网 络 的 层 内 耦合 强度 或 为 国定 ,或 为 可 变 。 
针对 多 层 变 耦合 星 型 网 络 的 同步 性 进行 了 研究 : 在 层 内 塌 合 强度 总 值 一 定 的 前 提 下 ， 令 网 络 中 部 分 耦合 强度 分 别 按 
有 照 正 比例 函数 、 反 比例 函数 、 指 数 函 数 、 对 数 函 数 以 及 三 角 函 数 等 规律 变化 ， 然 后 严格 推导 多 层 变 耦合 星 型 网 络 的 
超 拉 普 拉 斯 矩阵 及 其 特征 值 谱 ， 依 据 主 稳定 函数 法 ， 研 究 网 络 同步 能 力 。 数 值 仿真 结果 表明 : eit Ub duds 
度 发 生变 化 时 ， 其 余 es TM 可 能 出 现 不 同 区 间 内 网 络 同步 能 力 的 变 
化 趋势 不 同 的 情况 ; 另外 ， 同 步 域 有 界 时 网 络 同步 能 力 可 能 达到 极 大 的 值 。 综 合 分 析 可 以 得 出 台 EP 

的 变化 趋势 始终 与 较 小 的 耦合 强度 的 变化 趋势 相同 ; Rd E M 
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Research on synchronizability of multilayer star network with time-varying coupling strength. 
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Abstract: This paper studied the multilayer star network synchronization with variable coupling strength. First, from the 
premise that the sum of the intralayer coupling strength in the network is constant, part of the coupling strength are made to 
vary respectively according to the rule of different functions, then the Supra-Laplace matrix and its eigenvalue spectrum of 
multilayer star network are strictly derived. Finally, the network synchronizability is researched on the basis of the main 
stability method. The results show that when part of the coupling strength vary, the other coupling strength will vary in a 
opposite trend. Besides, the network synchronizability may be very large when the synchronous domain is bounded. The 
conclusion is that the trend of the synchronizability is always the same as the trend of the smaller coupling strength in the 
multilayer star network, and it is easier for the multilayer star network to be synchronous when the synchronous domain is 
bounded. 
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0 ”引言 链 路 相连 的 网 络 组 成 的 双 工 网 络 的 同步 性 ， 发 现 连 接 度 高 的 
D: 节点 更 有 利于 实现 双 工 网 络 同步 。 徐 明明 等 人 RP?”23 采 用 主 稳 
复杂 网 络 是 近年 来 较为 热门 的 研究 课题 "5]。 多 个 网 络 2 定 函数 法 深入 研究 了 两 层 网 络 的 同步 性 问题 ， 讨 论 了 不 同 耦 

间 可 能 存在 相互 联系 、 相 互 作 用 ， 从 而 形成 多 层 复 杂 网 络 ， 合 方式 和 网 络 结构 参数 对 网 络 同步 能 力 的 影响 ， 并 分 析 了 网 
更 能 精确 、 全 面 地 描述 现实 社会 中 复杂 系统 的 形态 、 作 用 和  ” 络 的 层 内 和 层 间 耦合 强度 分 别 较 弱 时 ， 如 何 使 网 络 更 好 的 实 
发 展 趋势 。 目 前 ， 多 层 网 络 已 经 在 运输 81、 社交 四、 经 济 四 。 现 同步 的 问题 。 孙 娟 等 人 中 进一步 研究 了 三 层 以 及 多 层 星 型 
等 诸多 方面 广 为 应 用 。 故 而 近年 来 对 复杂 网 络 的 研究 逐渐 向 网络 的 特征 值 谱 及 同步 性 ， 讨 论 了 网 络 耦 合 强度 、 层 数 、 节 
多 层 网 络 拓展 ， 已 经 在 拓扑 结构 0020、 动 态 网 络 模型 52 点 数 等 因素 对 网 络 同步 能 力 的 影响 。 
和 重 棒 性 52829 等 方面 取得 了 较 大 研究 成 果 。 其 实 ， 现 实 中 许多 复杂 问题 可 以 抽象 、 简 化 为 多 层 网 络 
同步 是 自然 界 中 一 种 常见 的 现象 ， 如 萤火虫 内 烁 同步 P、， 问题。 例如， 为 了 统筹 考虑 线 上 社交 网 络 和 线 下 社交 网 络 在 
掌声 同步 33、 路 由 器 同步 3] 等 。 在 现实 生活 中 ,同步 有 时 是 ”舆情 传 播 中 的 作用 ， 沈 乾 等 人 B9 首 次 提出 了 一 种 包含 “媒体 
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有 益 的 ， 有 时 却 是 有 害 的 9。 为 了 利用 有 益 的 同步 造福 人 类 层 一 线 上 层 一 线 下 层 ” 的 多 层 同 步 网 络 模型 ， 通 过 案例 仿真 
司 时 避免 有 害 同 步 ， 需 要 不 断 对 同步 问题 进行 深入 研究 。 随 分 析 线 上 网 络 与 线 下 网 络 之 间 的 相互 影响 及 其 对 与 论 传播 的 
着 多 层 网 络 研究 的 发 展 ， 对 多 层 网 络 同步 性 的 研究 已 取得 ] 作用 ; 李 延 龙 等 人 9 将 稀 玻 突 触 的 复杂 神经 网 络 抽象 为 多 层 
一 定 的 成 果 。Alberto 等 人 中 探讨 了 层 间 连接 模式 对 多 层 随 机 网络， 进而 研究 其 同步 中 耦合 参数 的 临界 特征 。 随 着 科学 和 
网 络 同步 性 能 的 影响 ， 发 现 层 间 连 接 模式 中 某 些 对 称 性 的 存 。 究 的 日 益 深入 ， 多 层 网 络 原理 将 会 应 用 到 越 来 越 多 的 实际 问 
在 往往 会 降低 整个 网 络 的 同步 能 力 。Li 等 人 ?9 研究 了 由 两 个 ” 题 当 中 去 。 
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众所周知 ,神经 网 络 中 神经 元 之 间 通 过 突 触 连接 。 最 近 ， 
Jungenitz 等 人 研究 发 现 : 在 神经 网 络 中 ， 重 复 的 电 刺 激 可 以 
通过 扩大 树 突 棘 来 加 强 突 触 ， 然 而 ， 这 个 过 程 中 琼 的 整体 大 
小 和 数量 没有 改变 ， 也 就 是 说 ， 当 刺激 加 强 了 一 组 突 触 ， 并 
且 它 们 的 树 突 棘 扩大 时 ,没有 同时 被 刺激 的 另 一 组 突 触 变 弱 ， 
它们 的 树 突 环 收 缩 B3。 类 似 的 情形 如 路 由 器 分 配 宽带 问题 ， 
当 茶 些 终 端 占用 较 多 宽带 时 ， 其 他 终端 占用 的 宽带 就 会 随 
减少 。 基 于 此 ， 本 文 研究 多 层 变 耦合 星 型 网 络 的 同步 性 : 在 
层 内 耦合 强度 总 值 一 定 的 前 提 下 ， 令 网 络 中 部 分 耦合 强度 分 
别 按 正 比例 函数 、 反 比例 函数 、 指 数 函数 、 对 数 函 数 以 及 余 
弱 函 数 的 规律 变化 《其 余 耦 合 强度 随 之 改变 ) ， 然 后 严格 推 
导 多 层 变 耦 合 星 型 网 络 的 超 拉 普 拉 斯 矩阵 及 其 特征 值 谱 ， 依 
据 主 稳定 函数 法 ， 研 究 网 络 同步 能 力 。 通 过 数值 仿真 结果 ， 
观察 不 同 变化 规律 下 网 络 同步 能 力 的 发 展 趋势 ， 进 而 分 析 总 


结 出 一 般 性 结论 。 


1 ”基础 理论 
M 层 网 络 中 第 久 层 第 i 个 节点 的 动力 学 方程 为 
d NS rium u KH, H K 
ao OD tat wi xD p x; Xx (1) 


Erg. 1<i<N，1<K<M ; xen" 是 第 KK 层 第 i 个 节点 的 状 
dE. fO: RoR 是 一 个 动力 学 函数 ;a" 为 第 KK 层 的 


层 内 耦合 强度 ; Ww e RV" HE KARRAR, 若 节 点 


i 与 节点 jGepDdHwi-ixe, MW; 反之 起 =0; HWE: 


N 


vi-- 3 wi ; Gd 7x0 RRE K EUST RRS REB ERAS 
FRAG ai" 为 第 天 层 与 第 瓦 层 的 第 个 节点 的 层 间 耦合 强 


度 ; Qf xD 表示 第 工人 个 节点 状态 变量 的 层 间 耦合 函数 。 
需要 指出 ， 网 络 属 于 哪 种 类 型 是 由 网 络 上 各 节点 的 具体 
动力 学 行为 f() 和 节点 间 有 具体 的 耦合 方式 决定 的 。 在 分 析 网 
络 同步 能 力 的 时 候 ， 可 以 将 网 络 结构 对 同步 的 影响 〈 即 外 耦 
BEEN IAE BS mg) 和 具体 动力 学 行为 及 节点 间 具 体 耦 合 
方式 对 同步 的 影响 分 离 出 来 。 也 就 是 说 ， 描 述 网 络 结构 同步 
能 力 可 以 只 考虑 外 耦合 矩阵 。 外 耦合 矩阵 可 以 采用 多 种 不 同 
的 表达 形式 ， 本 文选 用 拉 普 拉 斯 矩阵 。 
对 于 一 个 每 层 节点 数 都 为 N 的 M 层 星 型 网 络 模型 : a) 
Hi w*en"" 表示 第 必 层 的 层 内 连接 强度 矩阵 ， 用 si eR” X 
示 层 内 点 强度 矩阵 (对 角 和 矩阵 ) ， 则 第 天 层 的 层 内 拉 普 拉 斯 
和 矩阵 表示 为 "=s*-w* b) 假设 每 层 网 络 的 网 络 结构 都 相 
， 且 只 有 层 间 对 应 节点 相连 。 那 么 ， 用 Ww'e R” 表示 每 两 
层 网 络 的 层 间 连接 强度 矩阵 ，s'e RW 表示 每 两 层 网 络 的 层 
则 点 强度 和 矩阵 (对 角 和 矩阵 ) ， 则 每 两 层 网 络 的 层 间 拉 普 拉 斯 
WE L'-s'-w!'. 
整个 网 络 的 超 拉 普 拉 斯 矩阵 A=A:+A' ,其 中 A’ 为 层 内 
超 拉 普 拉 斯 矩阵 ，A' 为 层 间 超 拉 普 拉 斯 矩阵 ; 层 内 超 拉 普 拉 
斯 矩阵 为 


LT 


可 


zm 
ful 


L 0 0 
wu) tT |- Q) 
MOS 
层 间 超 拉 普 拉 斯 矩阵 为 
A =L QI (3) 


HEH: © Jy Kronecker 乘积 ,7 是 一 个 MxN 的 单位 矩阵 。 
求解 网 络 模型 的 超 拉 普 拉 斯 窍 阵 〈 外 耦合 矩阵 ) 的 特征 
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值 谱 。 根 据 主 稳定 函数 法 ， 模 型 的 同步 能 力 可 以 用 A 的 最 小 
非 零 特征 值 A2 和 最 大 特征 值 tex RAPSA, 当 同 步 区 域 
无 界 时 ，42 越 大 ， 网 络 同步 能 力 越 强 ， 当 同步 区 域 有 界 时 ， 
Ass A 的 值 越 小 ， 网 络 同步 能 力 越 强 。 


2 ”同步 性 研究 


2.1 网 络 结构 模型 

构建 M 层 变 耦 合 星 型 网 络 模型 ,每 层 网 络 共有 N 个 节点 
且 每 层 网 络 的 网 络 结构 和 耦合 方式 都 相同 。 图 1 (a) 为 网 络 
层 内 耦合 结构 图 ， 叶 子 节点 向 中 心 节点 单 向 耦合 ， 同 一 层 内 
叶子 节点 间 相 互 无 夺 合 ， 假 设 初始 时 刻 网 络 层 内 耦合 强度 为 
a 且 耦 合 强度 总 值 一 定 ， 为 (W-Da ， 某 时 刻 其 中 闫 个 耦合 强 
度 变 为 se ， 则 其 余 -mi ARa RR Jg 91, ， 为 方 


N —m-1 


El 


后 续 工作 ， 将 NCEU- 记 作 9。 图 1 (b) 为 网 络 层 间 耦 合 


结构 图 ， 每 两 层 网 络 的 对 应 节点 耦合 连接 ， 耦 合 强度 为 d, 
不 对 应 节点 间 无 耦合 。 


acc sa æ 
Oa . 
LAS Qa 0a L; 
Qa za; sa œ 
CS - 
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(a) 网 络 层 内 耦合 结构 


(a) Intra-layer coupling structure 


(b) 网 络 层 间 耦合 结构 ,红色 虚线 表示 中 心 节点 的 层 间 耦合 ,黑色 虚线 
表示 对 应 叶子 节点 的 层 间 耦合 
(b)Inter-layer coupling structure, the red dotted lines represent the 


coupling strength between corresponded hub nodes, the black dotted 
lines represent the coupling strength between the corresponded leaf nodes 
图 1 网 络 模型 耦合 结构 图 
(红色 节点 为 中 心 节 点 ， 黑 色 节点 为 叶子 节点 ) 
Fig. 1 Coupling structure of the network 


(the hub nodes are red and the leaf nodes are black) 
2.2 同步 能 力 
根据 前 述 理论 知识 ， 写 出 图 1 网 络 模型 的 层 内 拉 普 拉 斯 


JEREJ 
L*= 

(N -1)a*(M -1)d 0 m ds 25s 2 0 
一 EQ £a *(M -1)d : 
一 EQ ca *(M -1)d 
—0a Oa * (M -1)d 

: E 0 

E 0 0 86a«(M -1)d 


(4) 


ChinaXiv 合 作 期 刊 
录用 定稿 李晓霞 ， 等 : 多 层 变 耦合 星 形 网 络 同步 性 研究 第 37 卷 第 5 期 
民间 拉 普 拉 斯 矩阵 为 23 数值 仿真 
Ld da. 0 利用 MATLAB 进行 数值 仿真 , BLN 2200, M -20 , m-70, 
nolo : & ^' d=1， 仿 真 结果 如 下 : 
: "0 取 s=+ ， 对 网 络 层 内 耦合 强度 及 同步 能 力 进 行 仿真 如 图 
Vc ere ds oid 2 所 示 。 由 图 2(a) 可 以 看 出 ， 在 (0,2.84) 时 间 段 内 ， s 线性 增 
将 式 (4)(5) 代 入 到 式 (2)(3) 中 ， 分 别 得 到 层 内 和 层 间 超 拉 大 ， 而 4% 则 线性 减 小 ， 即 部 分 耦合 强度 增 大 ， 其 余 耦合 强度 
A Ay HE E, j tp E AS p 2 型 斯 矩阵 并 > "a a 
KREE, RENSMMETKUANONNERDIHUAUERORR A, ie pagg aca, Masa, raa, 
A B … B [ 
: B Ae B " 图 2(b)(c) RI AH A» 随时 间 增 大 ，7 随时 间 减 小 ， 即 网 络 同步 能 
其 特征 值 。 然 JA -|A- B[""|A« (M -DB 5d : 
特征 值 。 然 后 利用 公 或 |， :| SIAB MAOD 。 力 随时 间 增强 ， 在 (L.2.84) 时 间 段 内 ，a > ， 则 =9 ， 
B B Alum ela+Md v X — yy EM 
p &| 2(b)(c) 可 知 22 随时 间 减 小 ,随时 间 增 大 ， 
f: ER F IUE p 寺 征 信 分 别 为 0, ELUEA , £a-..0a ， Md...Md : ! 
推导 出 矩阵 的 所 有 特征 值 分 别 为 0， 867760 , 00-700 , Md- 有 
假设 Md> en, Md» Ga, 取 esc! MESA EAR GRIEF REJET CER 
则 最 小 非 零 特征 值 4.=minfsa94a} ， 最 大 特征 值 3 所 示 。 由 图 3(a) 可 以 看 出 ，s;: 非 线性 减 小 ， 而 9; 则 非 线性 
— —— EE rea MGe ， 此 可 得 增 大 ， 即 部 分 耦合 强度 减 小 ， "nc 在 (0.35,1) 
Wi ET À , 2265, 则 42-0», p EREMO 图 3(b H 
BUE IE RUE BERN F RE BUR i RU pu LES 
为 研究 网 络 时 变 耦 合 下 的 同步 能 力 ， 令 部 分 耦合 强度 分 4 和 :随时 间 增 大 ,> 随时 间 减 小 , 即 网 络 同步 能 力 随 时 间 增 强 ; 
m3: H8 TE 
BIHCRUIE LCS. RENE KAEN WAER E, m) MBRR, ncs Wase, rM, 
A aA E EE E E 7 
AROKISNCRMERUCOESUR eti eam ， e509 798 ， 图 eC sro Aa BENRA r 随时 间 增 大 ， 即 网 络 同步 能 
es=cost+1.5 , 力 随时 间 减 弱 。 
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(a) 层 内 耦合 强度 变化 曲线 (b) 同步 域 无界 时 ,X2 变化 曲线 (c) 同步 域 有 界 时 ,r 变化 曲线 
(a) Intralayer coupling strength (b) Graph of 42 when the synchronous (c)Graph of r when the synchronous 
region is unbounded region is bounded 
图 2 ”耦合 强度 按 正 比例 函数 变化 时 网 络 同步 能 力 


Fig.2 Graph of the synchronizability when the coupling strength vary according to the rule of proportional function 
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图 3 耦合 强度 按 反 比例 函数 变化 时 网 络 同步 能 


Fig.3 Graph of the synchronizability when the coupling strength vary according to the rule of inverse function 
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取 2-09 ， 对 网 络 层 内 耦合 强度 及 同步 能 力 进行 仿真 如 ”线性 减 小 , 即 部 分 耦合 强度 增 大 , 其 余 耦 合 强度 减 小 。 在 (1， 
图 4 所 示 。 由 图 4(a) 可 以 看 出 ，s; 非 线性 减 小， 而 0s 则 非 Ga+Md 
3 š sch ; 10) 时 间 段 内 ，s <2 ， 则 5s, r=. HR 5 œ) 
线性 增 大 , 即 部 分 耦合 强度 减 小 , 其 余 耦 合 强度 增 大 , 始终 都 i B ea 
4 随时 间 增 大 , r 随时 间 减 小 ， 即 网 络 同步 能 力 随 
fs. Maze, ro 271, HE 4C) (o) TAI Zo Bl Ce) 可 知 4: 随时 间 增 大 , 了 随时 间 减 小 Miedo ci 
" i ， 时 间 增 强 。 在 (10, ce ) 时 间 段 内 , e> 0s LI n m6, n E, 
间 减 小 ， 随 时 间 增 大 ， 即 网 络 同步 能 力 随时 间 减 弱 。 We Mo M MU oa 
取 esclogd ， 对 网 络 层 内 耦合 强度 及 同步 能 力 进行 仿真 由 图 5(b)(c) 可 知 A: 随时 间 减 小 ,> 随时 间 增 大 ， 即 网 络 同步 
如 图 5 所 示 。 由 图 5 Ca) 可 以 看 出 ， 4 线性 增 大 ,而 9: 则 ”能力 随 时 间 减 弱 。 
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(a) 层 内 耦合 强度 变化 曲线 (b) 同步 域 无 界 时 ，X; 变化 曲线 (c) 同步 域 有 界 时 ，7 变 化 曲线 
(a) Intralayer coupling strength (b) Graph of 22 when the synchronous (c) Graph of r when the synchronous 
region is unbounded region is bounded 
图 4 耦合 强度 按 指 数 函数 变化 时 网 络 同步 能 力 
Fig.4 Graph of the synchronizability when the coupling strength vary according to the rule of exponential function 
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(a) 层 内 耦合 强度 变化 曲线 (b) 同步 域 无 界 时 ，X2 变化 曲线 (c) 同步 域 有 界 时 ，r 变化 曲线 
(a) Intralayer coupling strength (b) The graph of 和 2 when the synchronous (c) Graph of r when the synchronous 
region is unbounded region is bounded 
图 5 耦合 强度 按 对 数 函数 变化 时 网 络 同步 能 力 


Fig.5 Graph of the synchronizability when the coupling strength vary according to the rule of logarithmic function 


Hy es cost L5 , TE(0, 4m) HT [8] Eit POSTE PT 288 E PE n ER EU 图 6(b)(c) RIA A» 随时 间 减 小 , ”随时 间 增 大 ， 即 网 络 同步 
同步 能 力 进行 仿真 如 图 6 所 示 。 以 一 个 周期 内 的 图 像 进行 分 。 能 力 随时 间 减 弱 。 
析 。 在 (0, ONARA, BA 6() 可 以 看 出 ，s; 减 小 而 0; 增 3 ”结束 语 
_ poeni 本 文 主要 研究 了 多 层 变 灯 合 星 型 网 络 的 同步 性 。 若 每 层 
As KB er 6s Roe. r7, 由 图 6(b)(9) 可 知 网 络 层 内 耦合 强度 总 值 一 定 ， 当 部 分 耦合 强度 分 别 按照 正比 
随时 间 增 大 ， 随时 间 减 小 ， 即 网 络 同步 能 力 随时 间 增 强 ， 例 函 数 、 反 比例 函数 、 指 数 函 数 、 对 数 函 数 以 及 三 角 函 数 的 
irony d 6 可 以 无 出。 避 先 减 小 后 增 大 。 规律 变化 时 ， 依 据 主 稳定 函数 法 ， 研 究 网 络 同步 能 力 的 变化 
f BRA, HIE GG)RTEUB IH, eOGHUNRGA aa, dca cles HINAMARCHRRUIULM. HASAR 
"EM ， 在 某 些 变化 规律 下 ， 可 能 出 现 看 
ea pi BERT DU IDEA AE. EXE 
Ju oe ouo. JUR eee iter "m 合 强度 交 蔡 变 小 的 情况 ， 导 致 不 同 区 间 内 网 络 同步 能 力 的 变 
图 6(bX(o) 可 知 妈 随 时 间 先 减 小 后 增 大 ，r 随时 间 先 增 大 后 减 ”化 趋势 不 同 ， 但 综合 分 析 发 现 ， 不 论 耦 合 强度 以 何 种 规律 变 
| ” 化， 不 论 同 步 域 有 界 还 是 无 界 ， 网 络 同步 能 力 的 变化 趋势 始 
小 , 即 网 络 同步 能 力 随时 间 先 减 弱 后 增强 ; E DE. ,2D) 时 间 上 l 
D HIPH FEREZUBUH MAREN ÆC DNAR 。 终 与 较 小 的 耦合 强度 的 变化 趋势 相同 ， 然 而 ， 同 步 域 无 界 时 
ORATOR, kreso, Maco, EM, 网络 的 同步 能 力 普遍 较 小 ,同步 域 有 界 时 网 络 同步 能 力 在 革 
Flo geo Ae T I Teen ba ^ 0 些 时 候 却 可 以 达到 极 大 的 数值 ， 由 此 可 知 ， 同 步 域 有 界 时 ， 
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网 络 可 以 更 快 地 实现 同步 。 
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(a) 层 内 耦合 强度 变化 曲线 
(a) Intralayer coupling strength 
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(b) 同步 域 无 界 时 ，X2 变化 曲线 
(b) Graph of 和 2 when the synchronous region is unbounded 
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(c) 同步 域 有 界 时 ，r 变化 曲线 
(c) Graph of r when the synchronous region is bounded 
图 6 耦合 强度 按 三 角 函 数 变化 时 网 络 同步 能 力 
Fig.6 Graph of the synchronizability when the coupling strength vary 


according to the rule of trigonometric function 
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